function [ProbabilidadHastaMaximo,DineroProbabildidad]=MontecarloxM(CT,combinacion,S,semanas) 
global numSimul
global numPredicciones
plotfigures = false;
%% 1 mes
dias=5*semanas;
 parfor simul=1:numSimul
    CT_1m(simul,:)=CT(simul,1:dias)-S;
 end

%plot Montecarlo Simulations
if plotfigures
 figure;
 hold on;
   for simul=1:1e4
        plot(CT_1m(simul,:), 'Color', [rand, rand, rand]);
   end
    title(strcat('Montecarlo xM: ', num2str(combinacion)));
    saveas(gcf, strcat('MontecarloxS_',num2str(combinacion)), 'jpeg'); 
end

   %% Probabilidad
   close all
   ajuste = 0;% Ajsutamos probabildiad si aumenta disminueye probabildiad y aumenta balance
   
   barras=100000;

 [nelements_1m,xcenters_1m]= hist(CT_1m(:,end),10000);
 [maximoDeelementos posicionmaximo]=max(nelements_1m);
 columnasMasprobabilidad = find(nelements_1m<=maximoDeelementos);
 
 %probabilidad de ganancias
 posicionCero=  min(find(xcenters_1m > 0));
 posicionmaximo = posicionCero+ajuste;
 ProbabilidadPerdidas = sum(nelements_1m(1:posicionCero))/numSimul;
 ProbabilidadHastaMaximo=1-(sum(nelements_1m(1:posicionmaximo))/numSimul);
 
 DineroProbabildidad= xcenters_1m(posicionmaximo);
 

 RangoInicial= min(xcenters_1m(columnasMasprobabilidad));
 RangoFinal= max(xcenters_1m(columnasMasprobabilidad));
 fprintf('Puede terminar(%i semanas) con un %i de probabilidad positivo \n',semanas,ProbabilidadHastaMaximo*100);

if plotfigures
        titulo = strcat('Puede terminar  con un ',num2str(ProbabilidadHastaMaximo*100));%'de probabilidad entre 'num2str(RangoInicial) 'y' num2str(RangoFinal)'\n'];
        titulo= strcat (titulo,'de probabilidad con ');
        titulo= strcat (titulo,num2str(DineroProbabildidad));
       hist(CT_1m(:,end),10000)
       title(titulo);
       saveas(gcf, strcat('Hist_MontecarloxM_',num2str(combinacion)), 'jpeg'); 
   end
 

end